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ABSTRACT 

Recent attempts to match shortwave albedo of snow for clear skies using approximate 
spectral solar fluxes and solutions of the radiative transfer equation for snow were unsuccessful 
until a separate surface reflection term was introduced. A separate consideration of specular re- 
flection from suiface snow grains has been objected to as being ad hoc. I discuss results based 
on a new parametcri/ation of shortwave radiation. Compared to the previous radiation models, 
the new model gives higher diffuse insolation and predicts higher albedos. The difference be- 
tween observed and predicted albedos is substantially reduced without invoking surface reflection. 


PPFCEDING PAGE BLANK NOT FILMED 


A Non ON SOI AR \ \ I VAI ION DI I’l NDI N('l Ol* 


a I'AR SKV SNOW Al HUH) 

The shoilwavo alhi'ilo of snow for dear skies was calciilatcii by H.itkstroin (1^72) using a 
solution of the railialive transfer et|uatii>n for a planc-patallel ‘atmosphere’ inoilel of snow. In 
this moilel ic..’ grains on the surface of snow scatter lailiation the same way as the grains ileep 
within. Hy showing gooJ agreement with observations (I iljeijuist, l*)5(i; Rusin, I *>(>4) he 
(piestioned atl Imc surface reflectivity introiluced by I iljet|uist (l‘>.^(>). Recently, WiscomlH* aiul 
Warren (1‘^HO) also have aigueil against this retlectivity. 

There were two shoitciunings in Haikstrorn’s calculation: (I) the ihffuse shortwave railiation, 
whose magnitude relative to ihiect raitiation increases as the solar elevation tlecieases, was igiuned, 
and (2) the albedos were not calculated in terms of snow grain si/e, but by adjusting a parameter 
in the radiative transfer ei|uatii>n, called the single scattering albedo. ('alculatii>ns of WiscomlH* 
aiul Warren (!*>S0) suffer the latter shortcoming. 

To remove these shortcomings, the albedos were calculated in terms of snow grain si/e using 
parameter i/ed models lor spectral solar tluxesand approximate M«lutions of the radiative transfer 
eijuation (I'houdhury, l‘>SI ;( houdhury and ('hang, l*>Sla, b). Tor a plane-parallel ‘alrrrosphere’ 
rrrotlel of snow, the calcrrlated albt*«los were lower than the obserxations. and the drffererree irrcreasi'd 
with ilecre.rsing si'lat elevatiorr (see ('hoirilhury, l*>8l ;('houdhury artd ('harrg, l‘>Slb). To match 
observatiorrs the surface retlectreity was reirrtrrMlitced. It was. however, cautioned that the failure 
to match observatiorrs withoirt a sc'parate surface retlectivity cv'itid be dire ti» inaccuracies in the 
rnoilels. This note gives albedo resirlts based on a new pararneteri/ation of incident Milar raitiation. 

The drffusi* railiation rnoilels useil previously were bawd irpiin the graphical rnetlunl sug- 
gested by Robinson (l‘>(»(>), which is to use an empirical ei|uation for the radiation on a surface 
with albedo 0.2.^, and correct it for the actual surface albeilo using his graphs for varied si'lar 



clrvulion&, iilinosphcric prccipiljbic water and turbidity. Ttic fact that snow is more reflective in 
the visible than in near infrared was not used in these models. One would expect that atmos- 
pheric backseat tering contributior depi'nds uptm spectral surface albedo. 

Sivkov (l‘>71) suggests that instead of calculating the diffuse radiation for varied surface al- 

K*dos, the total radiation be calcuKited from the erpiation 

4 4 3(1 - g) T 

g , . . 

4 4 3(1 - g) ( I - A) T 

where g is the effective asymmetry factor for atmospheric scattering phase functi<>n (Stdndev, 
l‘#75 and Mo et al., 1975), r is the sum of optical thicknesses for Rayleigh and aen>sol scatter- 
ings (leckner, 1978), A is diffuse allvdo of the surface (('houdhury and Hiang, 1979) and Qq is 
the total radiation vm a black surface i leckner, 1978 and Do/ier, 1980). Hie advantage of this 
luraineteri/ation over the Robinson’s is that spectral surface albedos can lx* inputed directly into 
the model. The ciH'ITicieni for (,\, can be shown to lx the atmospheric backscattering contribu- 
tion to total radiation in the delta-1 dilington approximation (loseph et al., 1976) of the radia- 
tive transfer equation. 

I'he overall accuracy of the radiation models may be assessed by comparing w’ith observed 

I 

iirsolations. llris comparison is shown in F'igure 1 for atmospheric parameters representative of 
Mirny and Maudheim (.Antarctic coastal stations); precipitable water 0.25 cm. Angstrom turbidity 
0.01, o/one 0.35cm, surface pressure 980mb and snow grain si^e 0.3mm. Hie model for direct 
radiation is iilentical to the previous imnlels (sec Choudhury and Oiang, 1981a), and the calcu- 
lated insolations are about 49( lower than the observations at high elevations. Hie present model, 
however, gives better results for diffuse insolations. 'Hie key factor in this improvement appears 
to be the inclusion of spectral snow' albedo in calculating the atmospheric backscattering contri- 
bution. It is to be noted that both direct and diffuse insolations now agree with the observations 
at low solar elevations where the observed albedos were found to iliffer most from the 
calculations. 
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I t^cw ;ilbciloi, together with those obtained using Choudhury and Chang's () 98 1 b) radiation 
model, are shown in Figure 2. The surface reflection is ignored. The present radiation model gives 
higher diffuse iso itions, and predicts higher albedos. The difference between the observed and pre- 
dicted albedos is reduced, but et low elevations the albedos are still underestimated by a few’percent. 
Whereas the calculated all>edos change very little at low elevations, the observed albedos change 
rapidly. Rusin (1964) also ob.ierved such variations. Approx bnate models should not be used to 
authenticate or discredit these obsen at ions. Rut, it appears unliVely that these variations are obser- 
vational errors associated with meas-uring low incident fluxes. The present radiation model agrees 
with the observed flux»*s, at least for low elevations. If surface gliid is a(l_ho£ and should be excluded 

from these calculations then further investigation of the accuracy of the dcIta-Fddington solution 

I 

(Joseph et al., 1976) for sm>w at low wlar elevations is needed to clarify these discrepancies. 
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Figure I. Comparison of calculated and observed insolations, (a) direct insolation, (b) diffuse insolation. The solid line using the 
present model, the dashed line using Choudhury and Oiang’s (1981b) radiation model. 






ALBEDO 



Fisurc Comp.inson ofol>s.*ncJ and calculated albedos. T1»e lines have the same 
nu^ning as in Figure Ih. Note the discrepancy >Mth observations at eloation 7.5 . 
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